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The Ray leigh- Taylor instability at the interface in an immiscible two-component Bose-Einstein 
condensate is investigated using the mean-field and Bogoliubov theories. Ray leigh- Taylor fingers 
are found to grow from the interface and mushroom patterns are formed. Quantized vortex rings 
and vortex lines are then generated around the mushrooms. The Ray leigh- Taylor instability and 
mushroom-pattern formation can be observed in a trapped system. 
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I. INTRODUCTION 



When a layer of a lighter fluid lies under that of a 
heavier fluid, the translation symmetry on the interface 
is spontaneously broken and the interface is modulated 
due to the Rayleigh- Taylor instability (RTI) 0, [2, @, Q. 
Waves on the interface then grow into complicated pat- 
terns with mushroom shapes yjllj. The RTI plays crucial 
roles in a variety of nonequilibrium phenomena, ranging 
from convection of water in a kettle to supernova explo- 
sions 

In the present paper, we investigate the RTI and ensu- 
ing dynamics in a phase-separated two-component Bose- 
Einstein condensate (BEC). Recently, there has been a 
growing interest in the interface properties of such BECs. 
For instance, the Kelvin-Helmholtz instability 0, 0, @], 
which occurs at the interface between two fluids with a 
relative velocity, has been observed in a 3 He superfluid 
system [9|, |l0j . The Kelvin-Helmholtz instability is also 
predicted in a two-component BEC of atomic gases . 
When a magnetic field is applied to a magnetic fluid 
(a colloidal suspension of fine magnetic particles), the 
surface is deformed by the Rosensweig instability [l2[ 
and grows into a pattern of crests. Such a surface phe- 
nomenon can be theoretically shown to occur also at the 
interface in a two-component BEC with a dipole-dipole 
interaction [l3j]. Analytical expressions of the interface 
tension in a two-component BEC have been derived in 
Refs. HE1. 



The present paper reveals that the RTI emerges at the 
interface between two immiscible BECs that are pushed 
toward each other by, e.g., a magnetic-field gradient. We 
first consider an ideal flat interface, and numerically show 
that the interface becomes deformed by the RTI to grow 
into the well-known mushroom pattern. The significant 
difference between this phenomenon and that in classical 
fluids is that the vortices under the caps of the mush- 
rooms are quantized. Thus, in three dimensions (3D), 
quantized vortex rings are generated around the mush- 
rooms. Bogoliubov analysis shows that the excitation 
spectrum of the interface modes closely resembles that 



for classical fluids. We also propose a realistic BEC sys- 
tem in a harmonic trap and show that the RTI can be 
observed experimentally. 

This paper is organized as follows. Section [HA] numer- 
ically formulates the problem and demonstrates the RTI 
and ensuing dynamics for an ideal system. Section iHBl 
gives the Bogoliubov spectrum of the interface modes. 
Section IIIII analyzes a trapped system. Section IIVI pro- 
vides conclusions to the study. 



II. RAYLEIGH-TAYLOR INSTABILITY AT AN 
IDEAL INTERFACE 

A. Mean-field dynamics 

The system considered here is a zero-temperature two- 
component BEC described by the Gross-Pitaevskii (GP) 
equations, 
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where x/jj, m^, and Vj are the macroscopic wave function, 
atomic mass, and external potential, respectively, for the 
jth (j = 1,2) component. The interaction parameters 
gjjf are given by 
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where a^y and m^y are the s-wave scattering length and 
reduced mass, respectively, between components j and 
j' . We assume that the interaction parameters satisfy 
the phase-separation condition, 
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FIG. 1: (Color) Density and phase profiles on the x-z plane. 
The initial state is a quasi- stationary state with a small ran- 
dom seed. The white arrows show the directions of the Stern- 
Gerlach force produced by the field gradient of B' = 20 
mG/cm. The red circle indicates the location of a topological 
defect under a mushroom cap. The density is normalized by 
4 x 10 14 cm" 3 . The field of view is 100 x 100 /mi. 



For concreteness, we employ the hyperfine states 
\F,rriF) = and |1,— 1) of a 8T Rb atom for com- 

ponents 1 and 2, respectively. According to the scatter- 
ing lengths reported in Ref. [17|, we have an = CL22 — 
100. 4<2b and avi = 101. 3<2b with ae being the Bohr 
radius, and then the condition for the phase separa- 
tion in Eq. (j3]) is satisfied. The spin-exchange dynam- 
ics |1,1),|1,— 1) — > 1 1,0), 1 1,0) can be suppressed by, 
e.g., the microwave- induced quadratic Zeeman effect [18|, 
which can lift the energy of the |1,0) state. When the 
hyperfine spins are parallel to the magnetic field £?, the 
magnetic- field gradient exerts forces ±/zbV|-B|/2 on the 
two components in opposite directions, where [ib is the 
Bohr magneton. 

In order to clearly demonstrate the RTI, we first con- 
sider a 3D system without a trapping potential. We 
prepare a quasi-stationary state with a field gradient 
B' = dB / dz > 0, which is uniform in the x-y direction. 
Components 1 and 2 are located in the z < and z > 
regions, respectively, and their fiat interface is located at 
the z = plane. The Stern- Gerlach force pushes the two 
components toward each other. We add small random 
seeds to the initial state, which numerically breaks the 
translation symmetry in the x-y direction and triggers 
the RTI. We assume periodic boundary conditions in the 
x and y directions. 

Figure [1] shows the time evolution of the density and 
phase profiles on a plane perpendicular to the initial flat 
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FIG. 2: (Color) (a) Isodensity surface of component 1 at 
t = 310 ms. The condition is the same as that in Fig. [1] (b) 
One of the mushroom shapes in (a). The red rings show the 
location of the topological defects. The red arrows indicate 
the directions of the atomic flow. 



interface. The two components are pushed toward each 
other (arrows in Fig. [I]), and the interface starts to modu- 
late due to the RTI (second row of Fig.[T]). Subsequently, 
the amplitude of the wave on the interface grows to form 
the mushroom shapes (third row of Fig. [1]). We can see 
that there are quantized vortices under the caps of the 
mushrooms (red circle in Fig.[T]). After that, the vortices 
enter into the mushroom patterns, giving rise to compli- 
cated dynamics (fourth row of Fig. [T]). 

Figure [2] (a) shows the isodensity surface of component 
1 at t = 310 ms. We can see that Ray leigh- Taylor fingers 
and mushroom patterns flourish at the interface. One of 
the mushroom shapes is magnified in Fig. [51(b), in which 
the topological defect under the cap of the mushroom is 
indicated by the red ring. The vortex rings are gener- 
ated by the upward flow of the atoms around the center 
and downward flow at the periphery of the cap of the 
mushroom shape. 



B. Bogoliubov analysis 

Before studying the Bogoliubov spectrum, we recall 
the dispersion relation for the RTI in classical fluids. Let 
us consider a situation in which inviscid incompressible 
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fluids produce a flat interface perpendicular to the di- 
rection of gravity, where the densities of the lower and 
upper fluids are p\ and p2 with pi < p2. From the lin- 
ear analysis, the interface mode has a dispersion relation 
of [| 
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where g is the gravitational constant and a is an 
interface-tension coefficient. If gravity is absent, uj is real 
for all k and proportional to k 3 ! 2 . In the presence of the 
gravitational force, there is always a range of k in which 
uj is imaginary, and hence the interface is dynamically 
unstable. The range of instability is given by 



< k < 



(P2 - Pi) 9 



2ttA: 



and the most unstable wave number fco is 
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We perform the Bogoliubov analysis by decomposing 
the wave function as 

Mr, t) = [fj(z) + 4> 3 {r)] e-w* (j = 1, 2), (7) 

where fj(z) is a quasi-stationary state with a flat in- 
terface around the z = plane and pj is the chemical 
potential. For the case of Fig. [TJ fi(z) and /^(z) are 
symmetric with respect to the interface on the x-y plane 
and pi = p2 = p because gu = #22- The small deviation 
(j)j{r) from the ground state is written as 

^(r) = u^ k (z)e^ kx -^ + vl k {z)e-^ kx -^\ (8) 

where the wave vector is assumed to be in the x direction 
without loss of generality. The mode functions Uj^(z) 
and Vj y k(z) satisfy the Bogoliubov-de Gennes equations, 
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where (3,3') = (1,2) and (2,1), m is the mass of 87 Rb, 
Vi = —p&B'z/2, V2 = pbB'z/2, and fj is assumed to be 
real. 

Figure [3] shows the Bogoliubov spectrum obtained by 
numerically diagonalizing Eq. (|9]). In Fig. [3l we plot 
only the lowest mode, which corresponds to the mode 
localized near the interface. The second lowest mode 
has much larger energy. When the field gradient B' is 
zero, the excitation energy is real for all wavelengths A 
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FIG. 3: (Color) Real part (solid lines) and imaginary part 
(dashed lines) of the Bogoliubov spectrum of an interface 
mode for a flat interface. The field gradient perpendicular 
to the interface is B' = (red), 5 (blue), and 15 mG/cm 
(green). The red line in the inset shows a logarithmic plot 
for B' — 0. For comparison, a black line with a slope of 3/2 
(u oc A -3 / 2 ) is shown in the inset. 



(red line in Fig. [3]). The logarithmic plot in the inset 
of Fig. [3] indicates that uj is proportional to /c 3//2 , which 
is in agreement with Eq. (|4]) with g = 0. In the pres- 
ence of the field gradient B' \ the Bogoliubov spectrum 
becomes imaginary for A -1 smaller than a critical value 
A~\ which is A" 1 ~ 2.1 x 10 4 m" 1 for B' = 5 mG/cm 
and A" 1 ~ 3.8 x 10 4 m" 1 for B' = 15 mG/cm. The long- 
wavelength modes are always unstable for B' 7^ 0, as in 
the RTI in classical fluids. 

The analytic expression of the interface tension in a 
phase-separated two-component BEC has been derived 
m Refs. [II [IE G3. For d\2 1 Q> — 1 <C 1, where a = an = 
<222, the interface tension a has the form, 



\/2w(ai2 - a), 
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where n is the atom density. Using the characteristic 
density 4 x 10 14 cm -3 for n and substituting Eq. ([TP]) 
into Eq. ([5|), in which (p 2 — p\)g is replaced by np^B 1 ', 
we obtain A" 1 ^ 2.1 x 10 4 m _1 for B' = 5 mG/cm 
and A" 1 ~ 3.6 x 10 4 m" 1 for B' = 15 mG/cm. These 
values of A" 1 are in good agreement with those in Fig. 03 
Using Eqs. (|6]) and ([TO)) , the most unstable wavelength 
is estimated to be c± 40 pm for B' = 20 mG/cm, which 
is in qualitative agreement with the wavelength of the 
interface modulation in Fig. [H 



III. DYNAMICS IN A HARMONIC TRAP 

We next consider a system confined in an axisymmetric 
harmonic potential Vtrap = m \^t z ( x2 + ^ 2 ) + ^yV 2 ]/ 2 - 
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FIG. 4: (Color) Column density f\^j\ 2 dy of the two- 
component BEC in a tight pancake-shaped trap with u xz — 
2tt x 100 Hz and u y — 2tt x 5 kHz. The number of atoms 
is N = 8.1 x 10 6 with an equal population in each compo- 
nent. The initial state is the ground state for B' — plus a 
small random seed. At t = 0, the field gradient is changed to 
B' — 1.3 G/cm. The white arrows show the directions of the 
Stern- Gerlach force for the two components. The unit of the 
density is (moo xz /h) 3 ^ 2 . The field of view is 64 x 64 /im. 

The radial and axial trap frequencies are uo xz = 2ir x 
100 Hz and u z = 2ir x 5 kHz, and the potential has a 
tight pancake shape. The initial state is the ground state 
of the GP equation for B' = 0, in which the interface 
is parallel to the x axis and components 1 and 2 (the 
hyperfine states |1,1) and |1,— 1) of 8T Rb) are localized 
in the z < and z > regions, respectively. A small 
random noise is added to the initial state to trigger the 
dynamical instability. At t = 0, the field gradient B' = 
1.3 G/cm is applied in the z direction. 

Figure 0] shows the time evolution of the column den- 
sity J \iljj\ 2 dy of each component, obtained by 3D simu- 
lation of the GP equation. At t ~ 80 ms, the interface 
starts to modulate with a wavelength ~ 8-10 /im due 
to the RTI. This wavelength is in qualitative agreement 
with ~ 8.7 /im estimated using Eqs. © and ([TO]) with 
the peak density n ~ 3.2 x 10 15 cm -3 . The modulation 
on the interface then develops into the mushroom pat- 



terns at t ~ 100 ms. After that, the system evolves in 
a complicated manner and eventually components 1 and 
2 are interchanged, localizing in the z > and z < 
regions, respectively. 

In Fig. [H one can see that both components invade 
around to the back of each other and the periphery of 
the condensate is disturbed. This is because the repulsive 
interaction between the two components is weak in the 
low-density periphery region, and the two components 
pass through each other. For a spherical trap, the pass- 
through phenomenon is more severe and a much larger 
number of atoms (> 10 9 ) is needed to clearly realize the 
RTI. If we use a square- well potential produced by, e.g., a 
flat-top beam [l9[ , the atomic density becomes more uni- 
form on the interface, which suppresses the pass-through, 
realizing an ideal RTI as discussed in Sec. Ill Al Increas- 
ing the inter-component repulsion [20[ using the Feshbach 
resonance can also suppress the pass-through. 

IV. CONCLUSIONS 

We have shown that the RTI and mushroom-pattern 
formation occur in a two-component phase-separated 
BEC, as in classical fluids. The significant difference be- 
tween the quantum RTI and the classical RTI is that 
quantized vortex lines and vortex rings are formed un- 
der the caps of the mushrooms. The Bogoliubov analysis 
showed that the excitation spectra of the interface modes 
are very similar to those of classical fluids. We proposed 
a possible experiment to observe the phenomena in a re- 
alistic trapped system. 

Various phenomena related to fluid instabilities may 
be reproduced in BECs with renewed interest. For ex- 
ample, splashing of drops and crown formation [2l| is 
considered to be related to instabilities which include the 
RTI. Droplet formation by the Plateau- Rayleigh instabil- 
ity [22[ in BECs is also an interesting future problem. 
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